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We derive upper and lower bounds on the fidelity susceptibility in terms of macro- 
scopic thermodynamical quantities, like susceptibilities and thermal average values. 
The quality of the bounds is checked by the exact expressions for a single spin in 
an external magnetic field. Their usefulness is illustrated by two examples of many- 
particle models which are exactly solved in the thermodynamic limit: the Dicke 
superradiance model and the single impurity Kondo model. It is shown that as far 
as divergent behavior is considered, the fidelity susceptibility and the thermodynamic 
susceptibility are equivalent for a large class of models exhibiting critical behavior. 
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I. INTRODUCTION 



Over the last decade there have been impressive theoretical advances concerning the 
concepts of entanglement and fidelity from quantum and information theory [l|,[2l,|3l, and 
their application in condensed matter physics, especially in the theory of critical phenomena 
and phase transitions, for a review see % IS] . These two concepts are closely related to each 
other. 

The entanglement measures the strength of non-local quantum correlations between par- 
titions of a compound system. So it is natural to expect that entanglement will be a reliable 
indicator of a quantum critical point driven by quantum fluctuations. To this end the main 
efforts have been focused on the different entanglement measures and their behavior in the 
vicinity of the critical point js-Sj. 



The fidelity jo-S] naturally appears in quantum mechanics as the absolute value of the 
overlap (Hilbert-space scalar product) of two quantum states corresponding to different 
values of the control parameters. The corresponding finite-temperature extension, defined 
as a functional of two density matrices, pi and p2, 

(1) 



1/2 1/2 



was introduced by Uhlmann 



-^(Pi,P2) = Tryp/ P2P1 
6] and called fidelit y by Jozsa 



by Jozsa, and used, e.g., by the authors of jlOMlSl]. is the square of J-'(pi,p2), however we 



Actually, the definition given 



shall adhere to the expression ffl, as most authors do. This functional has become an issue 
of extensive investigations 14l-ll8|. 



Closely related to the Uhlmann fidelity ([T]) is the Bures distance 19| . 

deipi, P2) = x/2-2^(pi,p2) (2) 

which is a measure of the statistical distance between two density matrices. The Bures 
distance has the important properties of being Riemannian and monotone metric on the 
space of density matrices. 

Being a measure of the similarity between quantum states, both pure or mixed, fidelity 
should decrease abruptly at a critical point, thus locating and characterizing the phase 
transition. Different finite-size scaling behaviors of the fidelity indicate different types of 
phase transitions 20|. The fidelity approach is basically a metric one 15| and has an 
advantage over the traditional Landau-Ginzburg theory, because it avoids possible difficulties 
in identifying the notions of order parameter, symmetry breaking, correlation length and so 
it is suitable for the study of different kinds of topological or Berezinskii-Kosterlitz-Thouless 
phase transitions. 

Due to the geometric meaning of fidelity, the problem of similarity (closeness) between 
states can be readily translated in the language of information geometry 3, Illl . The strat- 
egy here is to make an identification of Hilbert-space geometry with the information-space 
geometry. Note that fidelity depends on two density matrices, pi and p2, i.e. on the corre- 
sponding two points on the manifold of density matrices. On the other hand, being sensitive 
to the dissimilarity of states, it could be used to measure the loss of information encoded in 
quantum states. 

The above mentioned decrease in the fidelity J-'(pi,p2), when the state p2 approaches a 
quantum critical state pi, is associated with a divergence of the fidelity susceptibility Xi;'(pi) 



3 



22| which reflects the singularity of J-'{pi,p2) at that point. The fidelity susceptibility 
Xf{Pi), which is the main objects of this study, naturally arises as a leading-order term in 
the expansion of the fidelity for two infinitesimally close density matrices pi and p2 = pi + 6p. 
For simplicity, in our study we consider one-parameter family of Gibbs states 



p{h) = [Z{h)]-'exp[-f3H{h)], 



defined on the family of Hamiltonians of the form 



H{h) = T-hS, 



(3) 



(4) 



where the Hermitian operators T and S do not commute in the general case, /i is a real 
parameter, and Z{h) = Tr exp[— 13 H {h)] is the corresponding partition function. Note that 
the fidelity and fidelity susceptibility under consideration are defined with respect to the 
parameter h, including the important symmetry breaking case when the system undergoes a 
phase transition as h is varied. The fidelity susceptibility at the point /i = in the parameter 
space is defined as (see e.g. [l7|): 

-21n^(p(0),p(/i)) d'J^{p{0),p{h)) 



Xf(p(0)) := lini 



(5) 



h=0 



From ([2]) and ([5]) we obtain for the case of two infinitesimally close density matrices the 
following relation between the Bures distance and the fidelity susceptibility: 



dl{p{Q),p{h)) = XF{pmh^ + 0{h' 



(6) 



ity 



To avoid confusion, we warn the reader that for mixed states the definition of the fide 
susceptibility based on the Uhlmann fidelity ([T]), differs from the one derived in 
(see also [sD by extending the ground-state Green's function representation to nonzero 



18| 



23| . see also our Discussion. 



temperatures. This fact has been pointed out in 

The quantity ([5]) is more convenient for studying than the fidelity itself, because it depends 
on a single point h = an d does not depend on the difference in the parameters of the 
two quantum states 



17 



24| . Physically, it is a measure of the fluctuations of the driving 



term introduced in the Hamiltonian through the parameter h. In the case when pi and 
P2 commute, there are simple relations between the fldelity susceptibility ([5]) and thermal 
fluctuations. For examples, in the case of two states at inflnitesimally close temperatures the 



fldelity susceptibility is proportional to the speciflc heat 
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22| ; in the case of spin systems 
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described by Hamiltonian of the form (jl]), with S being a projection of the total 



spin 



;hen 



the fidehty susceptibihty is proportional to the initial magnetic susceptibility [16|, [171, 122| . In 
the non-commutative case, when the driving term does not commute with the Hamiltonian, 
such type of relation become more complicated. 

It is our goal here to derive lower and upper bounds on theoretical information measures 
like fidelity susceptibility (|5]) and, therefrom, the Bures distance ([6]), in terms of thermody- 
namic quantities like susceptibility and thermal average values of some special observables. 
At that we concentrate on the most interesting non-commutative case. 

The paper is organized as follows: in Section [III we introduce our basic notations and 
derive an expression for the spectral representation of the fidelity susceptibility in the general 
noncommutative case, see ( 1371) . which is convenient for the derivation of inequalities involving 
macroscopic quantities, like susceptibilities and thermal average values. In Sections IIIII and 
[IV] new upper, (HOl) . and lower, (H3|) . bounds on the fidelity susceptibility are derived. In 
Section |V] these bounds are checked against the exact expressions for the simplest case 
of a single quantum spin in an external magnetic field. We also consider the Kondo model 
(Section IVII) and the Dicke model (Section IVIII) as non-trivial examples for testing our upper 
and lower bounds, and drawing conclusions about the behavior of the fidelity susceptibility 
itself. Finally, in Section IVIIH we make some comments and compare our results to the 



inequalities derived in 



for the Green's function based definition of fidelity susceptibility. 



II. FIDELITY SUSCEPTIBILITY 



Here we consider the fidelity susceptibility ([5]) for the one-parameter family of Gibbs 
states (IH]) at the point h = 0. To this end we rewrite the density matrix ([3]) identically as 

p{h) := p(0) + V(0) + (l/2)/iV(0) + r,{h), (7) 



where 

rs{h) := p{h) - p(0) - V(0) - (l/2)/iV(0) (8) 
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is expected to be of the order 0{h?) as /i — > 0. Next we calculate directly the derivatives in 
Eq. ©: 



p'(0) = p(0)/3 
p"(0) = 



5(A)dA - (S)o 



5(A)dA- (5)o 



(9) 



+ ! d\S{\) [ dX'S{X') 



p1 p1 

- / dA / dA'5(A - \')Si\) - / dA(5(A)^)o 
Jo Jo Jo 



Here 



(10) 



(11) 



^(A) = e^'^"^e-^'^\ 

and (. . . )o denotes average value with the density matrix p(0). For the further consideration 
it is important to note that, in conformity with the normalization conditions Trp(/i) = 
Trp(O) = 1, we have 

Trp'(O) = Trp"(0) = 0. (12) 

Our aim is to use the expansion (JTj) to calculate the square root in the definition of fidelity 
([T]) with accuracy 0(/i^), as /i — 0, which is sufficient for obtaining the fidelity susceptibility, 
see Eq. ([5]). 

To introduce our notation and for reader's convenience we derive here the matrix rep- 
resentation for the fidelity susceptibility in some detail. By using a slight extension of the 
method of [131, we set 



pi/2(o)p(/i)pV2(o) = p(0) +X + Y + Z, 



(13) 



where X, Y and Z are operators proportional to h, h"^ and h^, respectively, to be defined 
below. By squaring Eq. f lT3|) . we obtain: 



p^/'^p{h)p^''^ =pX + Xp + X^ + pY + Yp + pZ + Zp + XY + YX + 0{h^ 



(14) 



Here and below, for brevity of notation we have set p(0) = p. Equating terms of the same 
order in h up to 0{h^) yields a set of three equations: 



VV(o)p'/' = pX + Xp, 

(l/2)/iV'/V"(0)p'/' =X' + pY + Yp, 
p^^^r3{h)p^/^ = pZ + Zp + XY + YX. 



(15) 
(16) 
(17) 
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To proceed with the calculations in the general case, when the operators T and S do not 
commute, we introduce a convenient spectral representation. To simplify the problem, we 
assume that the Hermitian operator T has a complete orthonormal set of eigenstates |n), 
T\n) = Tn\n), where n = 1,2, ... , with non-degenerate spectrum {T^}. In this basis the 
zero-field density matrix p(0) is diagonal too: 

{m\p\n) = pmSm,n, m,n = 1,2,.... (18) 

In terms of the matrix elements between eigenstates of the Hamiltonian T, Eq. (fT5|) reads 

hpimm\p'mn)pl/' = (p„ + Pn){m\X\n), (19) 

hence we obtain 

1/2 1/2 

{m\X\n) = /i(m|/(0)|n) ^" . (20) 

Pm + Pn 

Similarly, the second equation (fT6|l yields 

{l/2)h^pli^{m\p"{0)\n)pl/^ = (m|XV) + (p^ + p„)(m|y|n), (21) 

hence 

{m\Y\n) = _H^!^ + (l/2)h^{m\p"i0)\n) ^"' ^" . (22) 

Prn ~t~ Pn Pm ~l~ Pn 

Now we turn back to Eq. (fT3l) and take the trace of both sides: 



TrA/pV2(o)p(/i)pi/2(o) = 1 + TrX + TrF + TvZ. (23) 



13| . no relationship between TrX and TiY appears in our scheme. From 



In contrast to Ref. 
Eq. ( I20l) we obtain 

TrX = (V2)Trp'(0) = 0, (24) 

due to the explicit form of p'(0), see Qj. Next, taking into account Eq. f|T2|) and the 
expressions f l20p . (122 p for the matrix elements of X and Y, we calculate 

m,n ' m,n ' ' 

Thus, by substitution of the above results into fl23l) . we derive with accuracy 0{h'^) the 
following expressions for the fidelity ([1]), where pi = p(0) and p2 = p(/i), 

rMO).p(/,))^l-l/.-E """^''°""'^ (26) 
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and the squared infinitesimal Bures distance (jG]) 



4Mo),pw)^iA's: i<7'°'">| , (27) 

m,n 

Finally, the definition of (|5]) yields (see , [isl for a similar matrix representation of the 

quantum fidelity susceptibility) 

5; (28) 

Here matrix elements of p'(0) are given by 

(n|p'(0)|m) = ^!i:i^(n|^|m), m ^ n, (29) 



and 

(m|p'(0)|m) = pr^/3[{m\S\m) - {S)o\, (30) 

where 

(5')o:=^p„(n|5|n). (31) 

n 

With the aid of the above expressions, the fidelity susceptibility (128|) takes the form 



Pn - 

T -T 



H 2 



|(m|^|n)|2 1 



+ -/3^((5^'^)0o, (32) 



2 , p^ -|- p„ 4 

where 6S'^ = S**^ — {S'^)o, S'^ being the diagonal part of the operator S: 

■.= J2{MS\m)\m){m\. (33) 

m 

Hence, 

((5^')')o:=E/^™H^H'-(^)o- (34) 

TO 

Finally, by using the identity 

Pm + Pn = {Pn " Pm) COthX„„, (35) 

where 

Xmn ^ ^^^--^-\ (36) 

we rewrite Eq. fl521) in the form: 

xHp) = ^ E 'V''%'^iyj + ^/3^((^g'^)^)o. (37) 
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Here, the first term in the right-hand side describes the purely quantum contribution to the 
fidehty susceptibility, because it vanishes when the operators T and S comute, while the 
second term represents the "classical" contribution. 

Now we are ready to derive lower and upper bounds on the fidelity susceptibility by 
applying elementary inequalities for (xcothx)"^ to the summand in the above expression. 



III. UPPER BOUND 

The application of the inequality (xcothx)"^ < 1 to Eq. ( 137|) readily gives: 

/^^ Pn -Pm,/ I crl \|2 1 



(p)<^ E ^V^i<^i^i"^)i' + i^'«^^')')°- (38) 

m,n,m^n 

By comparing the right-hand side to the expression for the Bogoliubov-Duhamel inner prod- 



uct (see, e.g 



25 



26| and references therein) of the self-adjoint operator 6S with itself, 
{6S-6S)o:= J^'dr{6S{T)6S)o = l J] ^V^K^I^NI' + (('^5')')o, (39) 



we obtain an upper bound in the transparent form 

Xf{p)<^{6S;6S)o. (40) 

Note that the right-hand side of the above inequality is priportional to the initial thermo- 
dynamic susceptibility: 

{dS;6S)o = - — — \h=o.= jXn, (41) 

where f[H{h)] is the free energy density of the system described by the Hamiltonian and 
Xn is the susceptibihty with respect to the field h. 



IV. LOWER BOUND 

A lower bound follows by applying to the spectral representation for the fidelity suscep- 
tibility f l37|l the elementary inequality 

(xcothx)"^ > 1 - (l/3)x^ (42) 
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Then, from representation (137|) . we obtain the following lower bound for the fidelity suscep- 
tibility 

m,n,m,j^n 



Xf{p) > ^{SS;6S)o- — ^ {pn - pm)Xmn\{n\S\m)\' 



^iSS-SS)o-^{[[S,nS])o. (43) 



V. THE CASE OF A SINGLE SPIN IN MAGNETIC FIELD 

To test the quality of the derived upper and lower bounds, we consider the simplest 
case of a single spin in external magnetic field subject to a transverse-field perturbation. 
By choosing the magnetic field in the reference state along the 2;-axis, and the transverse 
perturbation along the x-axis, the perturbed Hamiltonian takes the form 

f3H{hi) = -hia^ - hsa', (44) 

where and are the Pauli matrices. This choice corresponds to (3H{hi) given by Eq. 
with 

PT = -ha', /3S = (45) 
The fidelity ([T]) now measures the dissimilarities between the density matrices 

Pih) = ^"P(^i^^ + ^3^^), (46) 
2 cosh ^yhj + hi 

and p(0). The average values taken with the density matrix p(0) will be denoted by the 
symbol (...)o- Let |1) and | — 1) be the eigenvectors of a', as well as of the operator 
/3T = —h^a': f3T\ ± 1) = T±i \ ± 1), where T±i = =F^3- Therefore, the density matrix p(0) 
is diagonal in that basis and its nonzero matrix elements are 

Next, for the first derivative (Q we have the expression 



P'(0) := ^/(huhs 



hi=0 



p [\%X)dX-{a^)o. (48) 
Jo 



Taking into account that (T^'| ± 1) = | =F 1), hence (cr^)o = 0, we obtain that 

(±l|p'(0)| ± 1) = 0, (±l|p'(0)| T 1) = (49) 
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Thus, for the fidehty susceptibihty we obtain from Eq. ( 128|) the following result 

^ _ \{-l\p'{Om\' + \Wm-l)\' _ tanh^ h .... 

^pyp) wn. — r~\ — ju — ■ 



2(p-i + Pi) 4/i 



Evidently, the same result follows from expression ( 132|) . where the only nonzero matrix 
elements are {-1\(3S\1) = {1\/3S\ - 1) = 1. 

Let us check now the upper bound ( HOl) . First, by setting in ( l39l) f3S = and (36 S = a^, 
we obtain 

P^SS;SS)„ = Ezl^^^.^gsm? + |i^|{l|^S| - 1)|^ = (51) 

il — i_l i_i — Ji tls 



Then, we note that f3S'^ = 0, hence 



tanh^ /icj tanh /i-? 



= ^)ir - ^ ^ 

which is equivalent to the initial inequality (xcothx)"^ < 1. 

Consider now the lower bound (1431) . An elementary calculation yields 

P'[[S,T],S] = -h,[[a\a%a']=4h,a\ (53) 

hence, 

/33([[5,T],5])o = 4/i3tanh/i3, (54) 

and 

tanh / 1 



This lower bound amounts to initial elementary inequality ( l42|) and it is always valid but 
makes sense for /13 < 3 only. 

VI. APPLICATION TO THE KONDO MODEL 

For a long time the Kondo model has been one of the challenging quantum many-body 
problems in condensed matter physics. Recent considerations have shown that solid state 
structures containing interacting spin systems (the Kondo model provides a specific example 
of such structures) are attractive candidates for quantum information processing 3|. 

In its simplest formulation, the effective Kondo Hamiltonian describes a single magnetic 
impurity spin interacting with a band of free electrons in a spatial domain A. The Hamilto- 
nian has the form: 

'H^ = Ho-J{Sin'' + S2ny + S3n'), (56) 
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where 

k,a 

is the Hamiltonian for the conduction electrons. The fermion operator c\ ^ {ck,a) creates 
(annihilates) a conduction electron with momentum k and spin a, J is the spin-exchange 
coupling between the magnetic impurity spin S = {Si, S2, S3), = s{s + 1), located at the 
origin r = and the conduction-electron spin densities, n(0) = {n^ ,n^) at r = 0. In this 
case the Hilbert space naturally takes the form of a tensor product Sj (g> C^''^^, where is 
the Hilbert space of the free electrons, and exp(— is of trace class and defines a Gibbs 
state Q- 

Let us place the system in a homogeneous external field h in the direction of the z-axis, 

n^(h) =71^- hS,. (57) 
This choice corresponds to H{h) given by (jlj) with 

T := n^, S := 5-3. (58) 
Due to the rotational invariance of the model Hamiltonian T, we have 

(^3)0 = 0, {Sl)o = {l/3)s{s + l), (59) 
where (...)o denotes average value with the density matrix 

p{0):- 



Tre-^^A ■ 

Now, our inequalities ( 140|) and (143|) yield the following bounds on the fidelity susceptibihty. 

^(^3; Ss)o - ^([[53, T],Ss])o < Xf{p) < ^(^3; ^3)0. (60) 

Here the double commutator and the Bogoliubov -Duhamel inner product in f l60|) can be 
estimated in terms of the model constants in ( |56l) with the aid of the remarkable inequalities 
obtained by Ropstorff {28 1: 

1) 0<([[^3,T],53])o<^Jtanh/3J, (61) 



and 



Xc ^ < PiSs; ^3)0 < Xc, (62) 
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where xc is the Curie susceptibihty (the initial susceptibihty of a free spin) 



and the dimensionless quantity 



(3s{s + l) 

Xc = ^ , (63) 



^ 2s(s + l) ^ ^ 



is introduced. 

Let us use the Ropstorff's inequahties in combination with our inequahties. From (160|) . 
(|6T|) and fl62|) we obtain the following bounds on the fidelity susceptibility 



4 . X 
> ^Xf{p) > Xc 



(65) 



e/3 3 

These bounds can be readily analyzed and the formula exhibits some interesting properties. 
Since by definition Xf{p) is a nonnegative function, the lowest meaningful bound in the 
left-hand side of f l65|) is zero. The upper bound gives the maximal possible value of Xf{p)- 
Note that f{x) := (1 — e~^)/x is a strictly monotonically decreasing function of x = [0, oo), 
with lim^^o /(x) = 1 and lim^-^oo f{x) = 0. As a result the equation 

1 „ 

-X = 

X 3 

has a unique solution x = x*. Therefore, the lower bound in fl65|l is nonnegative in the 
interval < f3e < x* . 
In the limit 

BJtanhBJ 
^ 2s(s + l) 

we obtain an asymptotic relation between the fidelity susceptibility and Curie magnetic 
susceptibility ( 163|) 

Xf(p) ^ -^Xc- (66) 

From the previous consideration it is known that the above relation becomes equality if 
the driving term S and the Hamiltonian T commute, irrespectively of the concrete model. 
Indeed, this is the case in the Kondo problem in the classical limit of infinite spin, s ^ oo, 
with /3 held fixed. Another possibility (|66l) to take place asymptotically is PJ — )■ 0. Note 
that, in the case oo > /3e > x*, where the lower positive bound follows from the definition 
of Xf{p), we obtain the inequalities 

< XFip) < jXc (67) 
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between the fidelity susceptibility and Curie magnetic susceptibility ( 163|) . 

In general, we see that [4/ I3)xf{p) can be smaller than the free spin susceptibility but 
exceeds some nonnegative value given by 



max 



1-e- 

7^ 



- -e/3, 
3 







which depends on the parameter /3e. 

If (3e is close to zero then (4//3)xf(p) is close to the free spin susceptibility. In other 
words, from the bounds obtained here it is seen that, when /3e — ?• 0, as well as when s — ?■ oo, 
the deviation of {4/ P)xf{p) from the single-spin Curie law tends to zero. 



VII. APPLICATION TO THE DICKE MODEL 



The Dicke model 29|] considers the interaction between N two-level atoms and the quan- 
tized single-mode electromagnetic field in a cavity A of volume V. Working at a constant 
density n = N/V, one can take the Hamiltonian in the form 

H§ = uja+a + eJ'^ + XN~^/'^{a + a+)J^. (68) 

Here a and a"*" denote the annihilation and creation operators for the cavity mode with 
frequency < u < oo which act on the one-mode Fock space J^b and satisfy the commutation 
relations [a"*", a"*"] = [a~,a"] = 0, [a"^,a~] = 1; = (1/2) ^^^af, where a" (a = x,y,z) 
are the Pauli matrices acting on a two dimensional complex Hilbert space C^, describe the 
ensemble of two-level atoms, the real e is the atomic level splitting and < A < oo is the 
coupling strength of the atom-field interaction. 

With the aid of different methods, it was rigorously proven that in the thermodynamic 
limit the Dicke model exhibits a second order phase transition from a normal to a superradi- 



ant phase 



3o|-32 



-32|. For values of the parameters A, u and e obeying the condition 4A^/a; < |e|, 
no phase transition occurs at any temperature, whereas for AX'^/u > \e\ there exists a finite 
critical temperature Tc = (|e|/2) tanh(|e|a;/4A^). At the point T = 0, 4A^/ci; = |e|, the 
model exhibits another, quantum phase transition of second order driven by the parameter 
A . Quite recently, a renewed interest to this issue appeared due to studies of quantum en- 
tanglament. The superradiant quantum phase transition at the critical point Ac = \/\e\uj/2 
manifests itself, among other phenomena, in the entanglement between the atomic ensemble 
and the field mode. 
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Let us introduce a source term 5* for the electromagnetic field in the form 

'H^(h) =V.^- hS, S := VN{a+ + a)/2, 

where h is a real field conjugate to the self-adjoint operator S. Then, the second derivative 
of the corresponding free energy density with respect to h is 

^-I!^!!:}l = -Pi5S;5S)n^,,,, (69) 

where SS := S— {S)'^D(^f^y Let us note that in the thermodynamic limit the second derivative 
diverges at the point (/i = 0, T = T^). 

In the case under consideration our inequalities P0|) and f H3|) yield 



^{^S- 6S)o - ^{[[S, n^], S])o < Xf{p) < ^{ss- SS)o, 

where the subscript "0" is used instead of the Hamiltonian H^. One readily calculates 

{[[S,n^],S])o = Nu, (70) 

and, by using the definition ( BTjl of the initial susceptibility per particle with respect to the 
driving field h, we obtain the bounds 

/3 xf(p) , (3 

4^^ ~ 1^ - ^ - 4^^- 

Using these inequalities one can estimate apart the classical contribution and the quantum 
contribution to the fidelity susceptibility. When the size of the system increases — ?■ oo 
at the critical temperature Tc, the thermodynamic susceptibility xn — ^ oo and for a fixed 
value of (5u/48 we have the asymptotical result 



N 4 



Xn- 



Such a relation is a feature of classical systems. The above relation should be exact if "H^ 
and S commute, which is not the case. However, one can state that due to the strong 
classical fluctuations there is a temperature region around with borders defined by the 
condition 

f^c'XN » 

where the influence of the quantum fluctuations on the relationship between xn and xf{p)/N 
is suppressed. 
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VIII. DISCUSSION 

The study of theoretical information properties of quantum models at the critical point 
is a hot topic of a primary interest in the last decade. First, it allows one to explore 
the quantum world as a resource in quantum information processing. Second, it sheds 
light on the critical phenomena from a different point of view, being an useful tool in 
the case when the standard Landau- Ginzburg approach based on the idea of symmetry 
breaking is hampered. As a result, the obtained so far knowledge opens an exciting possibility 
of the study in a unified framework of quantum information theory, critical phenomena 
paradigms, and novel quantum phases of matter. Our consideration is an additional step in 
this direction. 

Let us mention that important results on the role of critical fluctuations in systems 
with broken symmetries have been obtained in the past using famous inequalities due to 
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33l435l|. It should be 



Bogoliubov, Mermin- Wagner, Griffits, among others, see, e.g. 
stressed that a noteworthy property of this inequality approach is that the obtained results 
are exact and cannot be inferred from any perturbation theory. The benefits of having exact 
statements about strongly interacting systems are difficult to overestimate. A survey of 
some other inequalities, used to solve problems arising in the Approximating Hamiltonian 



Method, has been recently given in 



261 ] along with a number of references. 



Inequalities, as a tool for obtaining interesting exact statements, are also wide-spread 



and traditional in the information theory, see e.g. [ill, |36|, [3 

Our lower ( H3l) and upper ( HOl) bounds indicate that for the detection of phase transitions 
between mixed states the fidelity susceptibility per particle Xf/N is as efficient as the usual 
susceptibility x- This conclusion is in conformity with the commonly accepted view that 
quantum fluctuations are dominated by the thermal ones when Tj, > 0. However, one 
should keep in mind that our results were derived under rather restrictive conditions on the 
spectrum of the Hamiltonian. 



The same conclusion has been drawn in 18|, although for a different definition of the 
fidelity susceptibility. A number of authors, see e.g. sl, [l8|, start from the expansion of the 
zero-temperature definition of fidelity (at (3 = oo) 

\{Mh)\\Mm - 1 - ^/^V(oo), h^O, (71) 
where ijJo{h) is the ground state eigenfunction of Hamiltonian (jl]). In terms of the spectral 
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representation of the operator H{0) = T, this yields [2l|, |22|: 



Xf[oo} 



|(^.(0)|^|^o(0))P 



(72) 



This quantity is related to the imaginary-time correlation function 



GiS\T) = eir)[{SiT)S)o-{S) 



OJ) 



where 6 is the Heaviside step function, as follows 2l|, |22| : 



TGiS\T)dT. 



Its finite-temperature generalization is defined as 
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x3m 



/3/2 



TG{S\T)dT, 



(73) 



(74) 



(75) 



where G{S\t) is again given by fl73|) but the symbol (...)o now denotes the thermal average 
in the Gibbs ensemble with Hamiltonian H{0) = T. A subtle point in this definition is the 
upper integration limit (3/2, in lieu of /3, which was justified by path-integral arguments and 
a Quantum Monte Carlo method in [iSj. 



In Ref. 



23| it was found that definition ( 1751) of the fidelity susceptibility based 



on a definition of the fidelity itself, 



(76) 



which is different from Uhlmann's definition ([T]) that was accepted in the present study. 
Let us consider the bounds on x'piP) obtained in 18| : 

1 



(77) 



where (in our notation) 



ci.^(0,/3) = ^/3^((55T)o+E7^ 

n>m ^ ' 



\{m\S\n)\^ (1 



~2X„ 



Pn- 



l + e- 



-2X„ 



(78) 



Here (is^(0, /3)/i^ is the leading-order term in the expansion of the squared Bures distance 
between the states p(0) and p{h) as h 0. As expected from relation (|6]), expression ( ITHjl 
exactly equals our fidelity susceptibility Xf{p), as given by spectral representation (|32|1 . 
Thus, inequalities ( 1771) provide a new lower bound on Xf{p), 



Xf{p) > x3{(3) 



(79) 
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which is a kind of dynamical structural factor of the driving Hamiltonian S. An advantage 



of this bound is t 
series expansion 



lat it can be computed with the aid of the quantum Monte Carlo stochastic 



38|. From the spectral representation of 18| 



we conclude that the lower bound fl79l) follows from expression fl32|) and the elementary 
inequality 

At that, the resulting bound will be with the better coefficient 1/4 of the classical term 
instead of 1/8 as in ( ISUj) . 

We have established bounds on the fidelity susceptibility which are expressed in terms 
of thermodynamic quantities. An additional reasoning that stimulates this line of consid- 
eration of information-theoretic quantities, is that the experimental setup for measuring 
thermodynamic quantities is well developed. Thus, estimation of metric quantities with a 
thermodynamic-based experiment seems to be very appealing. We have shown that as far 
as divergent behavior in the thermodynamic limit is considered, the fidelity susceptibility 
Xf and the usual thermodynamic susceptibility x ^i-re equivalent for a large class of models 
exhibiting critical behavior. It remains for the future to study the effect of degeneracy, es- 
pecially macroscopic one, of the ground state on the upper and lower bounds for the fidelity 
susceptibility. 
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